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ABSTRACT 

Although WASP-14 b is one of the most massive and densest exoplanets on a tight 
and eccentric orbit, it has never been a target of photometric follow-up monitoring or 
dedicated observing campaigns. We report on new photometric transit observations 
of WASP-14 b obtained within the framework of ‘Transit Timing Variations @ Young 
Exoplanet Transit Initiative’ (TTV@YETI). We collected 19 light-curves of 13 indi¬ 
vidual transit events using six telescopes located in five observatories distributed in 
Europe and Asia. From light curve modelling, we determined the planetary, stellar, 
and geometrical properties of the system and found them in agreement with the values 
from the discovery paper. 

A test of the robustness of the transit times revealed that in case of a non-reproducible 
transit shape the uncertainties may be underestimated even with a wavelet-based error 
estimation methods. For the timing analysis we included two publicly available transit 
times from 2007 and 2009. The long observation period of seven years (2007-2013) 
allowed us to refine the transit ephemeris. We derived an orbital period 1.2 s longer 
and 10 times more precise than the one given in the discovery paper. We found no 
significant periodic signal in the timing-residuals and, hence, no evidence for TTV in 
the system. 

Key words: planets and satellites: individual: WASP-14 b - stars: individual: 
GSC 01482-00882 - planetary systems. 


* Based on observations collected at the Centro Astronmico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck Institut fiir Astronomie and the Institute de Astrof- 
sica de Andalucia (CSIC). Based on observations obtained with 
telescopes of the University Observatory Jena, which is operated 
by the Astrophysical Institute of the Friedrich-Schiller-University 
Jena. 

f E-mail:sraetz@cosmos.esa.int 


1 INTRODUCTION 

The star WASP-14 (V = 9.75 mag) w as found to have a 
transiting planet bv Ijoshi et al.l (l2009l ). The planet orbits 
the host star with a period of 2.24 d causing transits with a 
depth of 11 millimag (mmag) and a duration of 2.8 h. High- 
accuracy photometric and spectroscopic follow-up observa¬ 
tions allowed one to determine the mass of the planet Mp = 
7.34 ± 0.50 Mjup, the planet radius J?p = 1.28 ± 0.08 Rjup, 
and the orbital semimajor axis a = 0.036 ± 0.001 au. These 
findings show that WASP-14 b is one of the densest exoplan- 
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Table 1. Physical and orbital properties of the WASP-14 b sys¬ 
tem summarized from literature. 


Parameter 

Value 

Ref 

Epoch zero transit time Tq [d] 

2454963.93676 

[1] 


±0.00025 

[1] 

Orbital period P [d] 

2.2437704 ±0.0000028 

[1] 

Semimajor axis a [au] 

0.036 ±0.001 

[2] 

Inclination i [°] 

84.32 ±0.60 

[2] 

Eccentricity e 

0.087 ±0.002 

[31 

Argument of pericentre u) [°] 

-107.1±0.5 

[3] 

Mass star Ma [Mq] 

1.211 ±0.125 

[2] 

Radius star Ra [R^©] 

1.306 ±0.070 

[2] 

Effective temperature Teff [K] 

6475 ±100 

[2] 

Surface gravity star logpA 

4.29 ±0.04 

[2] 

Metallicity [-^j 

0.0 ±0.2 

[2] 

Mass planet Af^ [Mj^p] 

7.341 ±0.500 

[2] 

Radius planet R\^ [Rjup] 

1.281 ±0.079 

[2] 

Distance d [pc] 

160 ±20 

[3] 

Age [Gyr] 

~0.5-1.0 

[2] 

Spectral type 

F5V 

[2] 


References: 111 I.Tohn son et al] ll2009t). f2l|jo shi et all ll2009l ~). and 
[31 iBlecic eTTu ll2013h 


ets with an orbital period of less than three days. The stellar 
density, the effective temperature, the rotation rate, and the 
high lithium-abundance indicate a very young age between 
0.5 and 1.0 Gyr. All system parameters known from litera¬ 
ture are summarized in Table [T] 

A very in t erestin g feature of WA SP-14 b repor t ed by 
l.loshi et al.l l|2000l ~) and confirmed by iHusnoo"^ all (I 2 OIII I 
is its high orbital eccentricity (e = 0.087 ± 0.002) for its 
small orbital distance. At this distance, tidal interactions 
with the star are expected t o circularize the orbit of the 
transiting planet (Zah^llOT^). The nearly circular orbits of 
most close-orbiting exoplanets agree with tidal circulariza- 
tion time-scales sig nificantly shorter than the system age. 
iHusnoo et all (I 2 OII ) determined a circularization time-scale 
for WASP-14 b to be 5 X 10^ yr, which is significantly shorter 
than its age. The high eccentricity of WASP-14 b may thus 
indicate either a tidal circularization time-scale comparable 
or longer to the system age or the presence of an additional 
perturbing body in the system w hich may signifi cantly de¬ 
lay the process of circularization llMardling|l2007l l . It is also 
possible that the planet may have arrived on its current or- 
bit recently. 

Ijohnson et al.l (l2009h have found evidence for a spin-orbit- 
misalignment (A = —33.1° ± 7.4°) in the WASP-14 system 
by measuring the Rossiter-McLaughlin effect. 

WASP-14 b belongs to the class of highly irradiated hot 
Jupiters. Observations of the ther mal emission d uring three 
secondary eclipses with Spitzer bv iBIecic et al.l (I 2 OI 3 I 1 gave 
first indications about the atmospheric composition and 
thermal structure of the planet. The observations neither 
indicate a temperature inversion nor an effective heat ex¬ 
change between day and night side. The chemical composi¬ 
tion is consistent with the solar abundances. 

WASP-14 b is a very massive transiting planet. Over 50% 
of these massive hot Jupiters show a significant eccentricity 
while intermediate-mass planets tend to be in circular or¬ 
bits. Although this finding could be biased because the cir- 
cularization time-scales with planet mass dMardlind l2007l . 


circularization time-scale is longer for massive planets) it 
may suggest that the formation and evolution of these plan¬ 
ets differ from the scenarios for lower mass planets. There¬ 
fore, the observation of WASP-14 is particularly interesting 
to constrain theories of planet formation, migration, and 
planet-star interaction. 

Although the non-zero eccentricity and the spin-orbit- 
misalignment may indicate that there could be additional 
bodies in the system, WASP-14 b has never been a target of 
photometric follow-up observations or observing campaigns, 
dedicated to detec t and characterize s ignals of transit timing 
variatio n (TTV) . ISouth worthI ll2012ll re-analys ed the data 
sets of Ijoshi et al.l 1 20091 ') and Johnson et al.l (I 2 OO 9 I I and 
found a significantly different set of physical properties of 
WASP-14 compared to previous studies. For these reasons, 
WASP-14 b was selected as a target of our TTV campaign. 


2 OBSERVATIONS, DATA REDUCTION AND 
PHOTOMETRY 

First observations of WASP-14 were already carried out in 
2009 at the University Observatory Jena. Since 2011 Febru¬ 
ary WASP-14 b has been a target of our TTV campaign. 
Altogether we collected 19 light curves (LCs) of 13 individ¬ 
ual transit events. We used six telescopes (with apertures 
ranging from 0.25-2.2 m) located in five observatories dis¬ 
tributed in Europe and Asia. The participating observato¬ 
ries with their telescopes and instruments are summarized 
in Table (2] 

The photometric data were reduced by the standard proce¬ 
dures including subtraction of bias and dark and dividing 
by a sky flat-field. We calibrated the CCD images using the 
IRAI0 routines darkcombine, flatcombine, and ccdproc. 
Aperture photometry was performed with a modified version 
of the standard iraf routine phot. Differential magnitudes 
were calculated using the met hod of the optimize d artificial 
comparison star developed bv iBroeg et al.l (I 2 OO 5 I I. All avail¬ 
able field stars are combined to create an optimiszd artificial 
comparison star by taking a weighted average. Very faint or 
variable stars get a low weight while bright and constant 
stars dominate the artificial comparison star. The final LC 
is produced by the comparison of WASP-14 with this artifi¬ 
cial standard star. 

Ten different aperture radii were tested. The aperture that 
produced LCs with the smallest scatter (smallest root mean 
square, rms) was finally chosen. 

As a preparation for t he LC analysis we fitted the LCs 
with the JKTBBOP code (ISouthw orth et ^I2004 alf9i. which 
is ba sed on the ebop program ( Etzellll98ll : IPopper fc Etzell 
I 198 III . It allows us to remove photometric trends (caused 
by a colour difference between target and comparison star, 
differential extinction, changes in airmass during the obser¬ 
vations) by fitting polynomials up to fifth order simultane¬ 
ously to the modelling. Throughout this work the LCs were 
detrended by fitting a second-order polynomial. 

Finally the differential magnitudes were transformed into 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 2. Observatories and instruments which observed transits of WASP-14. 


Observatory 

Long. (E) 

n 

Lat. (N) 
(°) 

Elevation 

(m) 

Telescope Q) 
(m) 

Camera 

# Pixel 

Pixel scale 
("/pixel) 

University Observatory Jena 

11.48 

50.93 

370 

0.25 

CTK“ 

1024x1024 

2.23 





0.60 

STK*" 

2048 X 2048 

1.55 





0.25 

CTK-II 

1056x1027 

1.19 

Calar Alto Observatory 

357.45 

37.22 

2168 

2.20 

CAFOS“ 

2048 X 2048 

0.47 

Stara Lesna Observatory 

20.29 

49.15 

785 

0.50 

STIOXE 

2184x1472 

0.56 

Observatorio de Sierra Nevada 

356.62 

30.06 

2896 

1.50 

Vers Array: 



(OSN) 





2048B 

2048 X 2048 

0.23 

TUBITAK National Observatory 

30.34 

36.82 

2485 

1.00 

SI 1100 Cryo 

4096x4097 

0.31 


“iMugraueJ ll2n09h . ^Mugrauer fc Bertholdl ll20inh . iMeisenheimerl lll994h 


Table 3. Summary of the WASP-14 b observations in the period from 2009 April to 2013 April: A^exp ^ number of exposures, Texp - 
exposure times, P — median number of exposures per minute, pnr - photometric noise rate. 


Date 

Epoch“ 

Observatory^ 

Filter 

-Vexp 

Texp (s) 

rms (mmag) 

F 

pnr 

2009 Apr. 01 

205 

Jena-CTK 

/ 

324 

60, 50, 40 

5.16 

0.87 

5.54 

2009 Apr. 19 

213 

Jena-CTK 

/ 

320 

30 

5.16 

1.02 

5.11 

2009 May 07 

221 

Jena-CTK 

R 

320 

25, 20 

10.03 

1.22 

9.08 

2011 Feb. 12 

509 

OSN 

R 

1215 

10 

3.56 

0.71“ 

4.23 

2011 Mar. 02 

517 

Jena-STK* 

R 

395 

30 

3.27 

1.39 

2.77 



Jena-CTK-II* 

R 

366 

40 

5.00 

1.39 

4.24 



Calar Alto* 

B 

144 

30 

6.07 

1.09 

5.83 

2011 Mar. 11 

521 

Jena-STK 

R 

330 

45 

2.46 

1.10 

2.34 

2011 Mar. 20 

525 

Jena-STK 

R 

360 

30 

1.87 

1.39 

1.59 



Jena-CTK-II 

V 

535 

30, 25 

4.65 

2.17 

3.16 



Calar Alto 

B 

322 

30 

10.36 

1.10 

9.87 

2011 Mar. 29 

529 

Jena-STK 

R 

401 

30 

1.25 

1.39 

1.06 



Stara Lesna 

R 

696 

30 

3.19 

3.02 

1.84 

2011 Apr. 07 

533 

OSN 

R 

570 

30 

4.10 

0.34“ 

7.08 

2011 Apr. 16 

537 

Calar Alto 

V 

404 

25 

9.99 

1.20 

9.13 

2012 Feb. 06 

669 

Jena-STK* 

R 

310 

30 

4.21 

1.39 

3.57 



Jena-CTK-II* 

V 

208 

60 

6.58 

0.95 

6.75 

2012 Feb. 24 

677 

Calar Alto 

V 

350 

30 

0.90 

0.39“ 

1.44 

2013 Apr. 30 

869 

TUBITAK 

R 

294 

30 

3.34 

0.75 

3.85 


‘partial transit or gaps i n the LCs. _ 

“calculated using the ephemeris in ijoshi et al.l ||2009| ~). 

*’for a description see Table[2] For the University Observatory Jena also the used camera is listed to avoid confusion. 

“for the final binned LC. 


fluxes and divided by the average out-of-transit value in or¬ 
der to normalize the LCs. All 19 LCs of WASP-14 are shown 
at the end of the paper. To give an estimate of the varying 
quality of the o bserved LCs, we com puted the photometric 
noise rate fonr: Ihlilton et aI]|2011^ . The pnr is calculated 
by dividing the rms of each LC, which is a result of the 
LC fitting with jktebop, by the square root of the median 
number of exposures per minute (L), 


A summary of all observations is given in Table [S] 


2.1 University Observatory Jena photometry 

Most observations were carried out at the University Obser¬ 
vatory Jena which is located close to the village Grofischwab- 
hausen, 10 km west of the city of Jena. There we have 
three telescopes (90, 25, and 20 cm) on the same mount 


each equipped with an optical CCD camera. For our tran¬ 
sit observations we used the 90 cm Schmidt telescope (60 cm 
free aperture in Schmidt-mo de) with the CCD-camera STK 
dMugrauer fc Bertholdl[2010l l and the 25 cm Casse grain aux¬ 
iliary telescope with the CCD-camera CTK dMugrauer! 
1200911 . In 2010 August, the CTK was replaced with a new 
CCD camera called CTK-II {Cassegrain Teleskop Kamera- 
II). The properties of the camera are given in Table ID 
Between 2009 March and 2012 March, we observed 11 LCs 
of eight individual transits of WASP-14 at the University 
Observatory Jena. Six of these eight transit are fully cov¬ 
ered. The remaining two transits show gaps in the data due 
to passing clouds. For the transit from 2012 February 6 only 
a part of the flat bottom is missing, so that the shape of 
the transit is still identifiable. The situation is different for 
the transit from 2011 March 2. In this case the whole egress 
was lost due to clouds that makes it difficult to determine a 
precise transit time. 

Observations were performed in Johnson V, Cousins R or 
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Table 4. Camera facts of the CTK-II 


Parameter 

Value 

Detector: 

E2V CCD47-10 

Pixel: 

1056 x 1027 

Pixelscale: 

(1.1956 ± 0.0001) arcsecpixeU^ 

Field of view: 

21.0 arcmin x 20.4 arcmin 

Filter: 

Bessell V, R, /, Gunn 2 : 

Focus: 

Cassegrain 


I filters. The optics were defocused slightly for the observa¬ 
tions with the STK. The exposure times varied from 20 to 
60 s with the 25 cm Cassegrain telescope and from 30 to 45 s 
with the 90/60 cm Schmidt telescope, depending on weather 
conditions, airmass, and telescope defocusing. 

2.2 TTV@YETI photometry 

The observation of a large number of individual transit of 
a planet in front of its host star by a single observatory is 
difficult. Due to the weather conditions in Central Europe 
and the orbital period of the transiting planet, it is almost 
impossible to observe consecutive transits. 

In 2009, we launched an international observation campaign 
which is dedicated to the detection and characterization 
of TTVs for carefully selected transiting planets. The pro¬ 
gramme is realized by the observation with globally dis¬ 
tributed telescopes at different longitudes and is based on 
the cooperation in the fr amework of the ‘Young Exoplanet 
Transit Initiative’ fYETI. lNeuhauser et al.ll201ll i. A descrip¬ 
tion of ‘TTV® YETI’ is available at the website of the 
projeciQ. 

2.2.1 Calar Alto 2.2 m telescope 

Between 2011 February and April, we awarded 2.5 nights 
(5x0.5 nights, project Fll-2.2-007) with the Calar Alto 
Faint Object Spectrograph (CAFOS) at the Calar Alto Ob¬ 
servatory. Due to extremely bad weather in winter 2011 two 
out of five transits were lost completely and one was only 
observed partially. The remaining two events could be ob¬ 
served but also under bad conditions like thin clouds, fog, 
and full moon. Therefore, the data have a relatively low qual¬ 
ity which is not sufficient to measure precise transit times. 
One additional LC was observed as back-up of project F12- 
2.2-009 on 2012 February 24. 

For the observations we used CAFOS in imaging mode 
and 2x2 binning. We windowed the field of view to 
7.9 arcmin x5.2 arcmin to shorten the read-out time. The 
observed field included WASP-14 and three suitable com¬ 
parison stars. Because WASP-14 is a relatively bright star 
the observations were carried out in Johnson B or V band. 
It turned out that the best LCs could be observed in the V 
band. To minimize random and flat-fielding errors, we defo¬ 
cused the telescope significantly. The telescope was autogu- 
ided during the observations. Depending on the atmosphere 
transparency and telescope defocusing, we used exposure 

^ http://www.home.umk.pl/ gmac/TTV 


times of 25 or 30 s. For the observations from 2012 February 
24, we combined three measurements into one data point to 
obtain a binned light curve. 


2.2.2 Observatorio de Sierra Nevada 

Two complete transit LCs (2011 February 12 and 2011 
April 7) were obtained at the Observatorio de Sierra 
Nevada using the 1.5 m reflector. With a VersArray:2048B 
CCD camera with 2048 x 2048 pixels and a pixel scale 
of 0.23 arcsec pixel”^, we could observe a field of view of 
7.85 arcmin x 7.85 arcmin. The exposure times of the R- 
band observations were chosen between 10 and 30 s. While 
the first LC were obtained under good weather conditions, 
the second transit suffered from passing clouds in the first 
half of the observation. To reduce the scatter, we binned the 
LCs by averaging every five data points. 


2.2.3 TUBITAK National Observatory 

One transit that occurred in 2013 April 30 was observed 
using the Spectral Instruments SlllOO series 4096x4096 
CCD camera mounted on the 1.0 m Telescope (TlOO) at 
TUBITAK National Observatory (TUG) in Turkey. The im¬ 
ages were acquired in R band with an exposure of 15 - 30 s 
under good weather conditions. 


2 . 2.4 Star a Lesnd Observatory 

One additional LC of WASP-14 was observed at the Stara 
Lesna Observatory in Slovakia. The observation in the night 
of 2011 March 29 was performed with a 50 cm Newtonian 
telescope equipped with an SBIG STIO MXE CCD camera. 
The CCD-chip consists of 2148 x 1472 6.8 pm pixels and 
has a pixel scale of 0.56 arcsec pixel”^ corresponding to the 
field of view of about 24x16 arcmin.The LC was obtained 
in Cousins R band with an exposure time of 20 s. 

Unlike all other observations used in this study, the data re¬ 
duction and photometry for this transit was already carried 
out at the Stara Lesna Observatory. The standard correction 
procedure (bias, dark and fiat-field correction) and subse¬ 
quently aperture photometry was performed by c-munipack 
packag^. Since the c- munipack pac kage is also based on 
the DAOPHOT program (IStetson|ll987l i like iraf, the results 
of both routines are comparable. To generate an artificial 
comparison star, at least 20 - 30 % of stars with the low¬ 
est LC scatter were selected iteratively from the field stars 
brighter than 2.5-3 mag below the saturation level. To mea¬ 
sure instrumental magnitudes, various aperture radii were 
used. The aperture which was found to produce LC with 
the smallest scatter was applied for generation of final LC. 
Due to problems with the flat-field and because the obser¬ 
vations were acquired without autoguiding the LC contains 
correlated noise. 
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Table 5. System parameters resulting from the simultaneous 
wavelet-based red noise MCMC analysis of the four best qual¬ 
ity LCs. 


Parameter Value 


Period 

2.2437704* 

Inclination 

85.3 

a/RA 

o.yo -0.32 

Rb/RA 

0 0965 

u.uyoo _o.oo27 

Linear LD (R band) 

0 39 +0.29 
u.oa _o.24 

Quad LD (R band) 

n 1 Q +0.38 
-0.45 

Linear LD (V band) 

Q CO +0.28 

-0.30 

Quad LD (V band) 

0 12 +0-41 
-0.40 

Eccentricity 

0.087* 


‘Value fixed in MCMC analysis. 


3 LIGHT-CURVE ANALYSIS 

To refine the parameters of the system and to determine pre¬ 
cise mid-transit times, it is necessary to model the individ¬ 
ual LCs. Therefore, we used the Transit Analysis Packag^fl 
(tap v2.1: lGazak et aDl2012l i. which employs Markov Chain 
Mo nte Carlo (MCMC) te chniques to fit transit LCs using 
the iMandel fc Agoll (|2002h model. To calculate the model, 
TAP uses the fast exoplanetary fi t ting c ode exofast which 
was developed bv lEastman et aP (l2013l l. tap has some ma¬ 
jor adva ntages. First, it incorp orates the wavelet-based tech¬ 
nique of ICarter fc Winra (l2009l i which allows one to estimate 
more robust parameter uncertainties than classic meth¬ 
ods because it parameterizes uncorrelated as well as time- 
correlated noise. Secondly, it is able to simultaneously anal¬ 
yse multiple transits observed in different conditions (instru¬ 
ment, filter, weather, etc). Finally, the tap code employs the 
quadratic limb darkening (LD) law which is a better choice 
to represent the obser vations than a sim ple linear law as 
shown, for example, in iRaetz et al.l ll2014ll . The theoretical 
LD coefficients that were used as initial values for the fit¬ 
ting were bilinearly interpolated (in ef fective tempe rature 
and surface gravity) from the tables bv IClare 3 (l2000li using 
the stellar parameters in Table [T] 

The LC analysi s was carrie d out as explained in 
iMacieiewski et ^ (l2013l . a) and iMacieiewski et ^ (l2013l . 
b). We selected the best-quality LCs sorted by their pur 
for a simultaneous fit using 10 MCMC chains with 10® steps 
each. The orbital inclination i , the semimajor-axis scaled by 
stellar radius , and the planetary to stellar radii ratio ^ 
were linked together for all LCs. We also connected the LD 
coefficients u and v but only for the LCs that were observed 
in the same filter. The orbital period P, the eccentricity e 
and the argument of periastron cj were kept fixed while the 
mid-transit times Tc, the airmass slopes, and the flux offsets 
were allowed to vary separately. The selection of the best 
LCs was done iteratively. The fitting procedure started with 
a few LCs with small pur and was repeated several times 
after adding more LCs with higher pnr. The final selection 
consists of the four best LCs with a pnr < 2 mmag, three 
LCs in R and one in the V band. Including transits with 

® http://c-munipack.sourceforge.net/ 

^ http://ifa.hawaii.edu/users/zgazak/If A/tap. html 



Figure 1. Phase-folded LC of the four best LCs with a pnr < 2 
mmag that were used to create the template shown as solid black 
line. The dark grey points show the same fluxes binned in phase, 
with a bin size of 0.005 d. 

a higher pnr degraded the quality of the fit. A phase-folded 
LC of these four transits including the best-fitting model 
that correspond to the resulting system parameters given in 
Table [5] are shown in Fig. [T] 


4 PHYSICAL PROPERTIES 

The results given in Table [5] allow us to calculate stellar, 
planetary, and geometrical parameters. The mean stellar 
density pA can be derived directly from the parameters ob¬ 
tained from the light-curve modelling using 

i|w in nl l20inll . where G is the gravitational constant. To de¬ 
termine the stellar mass we used the from Table [T] and 
the stellar density to plot WASP-14 into a modified version 
of the Hertzsprung-Russel diagra m (HRD), together with 
PARSEC isochrones (version 1.2S: lBressan et al.ll2012l l. The 
result for = 0.0 is shown in Fig. [2l WASP-14 is in an 
area of the HRD with overlapping isochrones of very young 
20 Myr) and young (~ 1 Gyr) ages. Since Ijoshi et al.l 
200ilh estimated an age range of 0.5-1.0 Gyr based on 
lithium abundance, the rotation v e locity and a comparison 
with the models of iFortnev et al.l (l2007tl . we excluded the 
very young ages. Taken also the uncertainty in the metal- 
licity into account the stellar evolutionary models yielded 
a stellar mass of Ma = 1.30 ± 0.06 Mpi which is consistent 
with the mass given in ijoshi et al.l (l2009ll . The improved pe¬ 
riod (see Section [5]), the stellar mass, the orbital Inclination 
i, the eccentricity e, and the amplitu d e of t he star’s radial 
velocity Ka taken from iBlecic et al.l ll2013ll allowed us to 
redetermine the planetary mass Mb. To calculate the semi¬ 
major axis a, we inserted the orbital period and the masses 
of star and planet into Kepler’s third law. By resolving a/ Ra 
and Rb /Ra using the already determined value for a we de¬ 
duce values for Ra and Rh- The planetary radius as well as 
the mass were used to calculate the planetary density pb- 
The impact parameter h which depends on the eccentricity 
e were calculated using 
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Figure 2. Position of WASP-14 in the plane. The 

PARSEC isochrones of solar metallicity for log(age) = 7.24-7.30 
with steps of 0.01 and log(age) = 8.95 - 9.40 with steps of 0.05 
for the very young age and the young age, respectively, are also 
shown. 


a cos i 1 — 

Ra 1 + e cos u 


(3) 


The surface gravities of star and planet qa and were cal¬ 
culated using the formulae 


9b = 

/27r/ 

(1-eY'^ 

a^KA 

(4) 

[p J 

' sinz 

Rl 

and 





9 A = 

f2TV' 


a^Ki, 

(5) 

[p , 

1 sin i 

Rl 

with 





Ab = 


2naMAami 


(6) 

(Ma ± Mb)P^ 

— 


l|Southworthl 120091 1. where Ka and Kh are the ampli¬ 
tude of the star’s and planet’s radial velocity, respec¬ 
tively. In addition, we calculated the equilibrium temper¬ 
ature of the planet Teq (assuming a Bond albedo = 0 and 
only li ttle energy redistribution across the surface of the 
pl anet; Hansen fc Barmanll2007l l and the Safronov number 
0 llSafronovI 1972l l. the square of the ratio of escape velocity 
of the planet v^sc and orbital velocity Uoru: 


1 / Uesc \ ^ _ a Mb 

2 VUorb/ Rh Ma 


(7) 


The results of the calculations are summarized in Table El 
Our r esults are consistent with the values in Ijoshi et all 
(l2009ll but have larger error bars. This is not surprising tak¬ 
ing into account the qualit y of our ground-ba sed LCs. The 
physical properties given in ISouthworthI (1201211 . however, are 
significantly different from our findings. The reason for this 
discrepancy lies in the interpr etation of the results of the 
LC analysis. ISouthworthI (l2012l l modelled two available LCs 
and obtained a set of solutions (for different LD laws and 
for fitting either one or two LD coefficients) which was av¬ 
eraged to a final solution. Since both LCs yielded different 
results, the final derived averaged properties deviate from 
the already published values. 


Table 6. Physical properties of the WASP-14 sys t em d erived 
fro m LC modell i ng. V alues derived by l.loshi et al.l ll2009l . J09) 
and ISouthworthI 120121 . S12) are given for comparison. 


Parameter 

This work 

J09 

S12 

Planetary parameters 

-^b [R-Jup] 

1 240 + ®-^^® 

i.z^u m Q JQ 3 

1 OQ"! _L 0.075 
1.251 ± 0.082 

1.633 ±0.092 

Mb [Mjup] 

7 CQ _|_ 0.24 
i.oy diQ 23 

7.34 ±0.50 

7.90 ±0.46 

Pb [pjup] 

3.73 ± Qigg 

o CA _|_ 0.64 

5.0U It 0 50 

1.69 ±0.25 

logSb 

4.090 ±0:080 

4 ni f) _{_ 0.049 

4i.UlU It 0.042 

3.866 ±0.042 

^eq [K] 

1872 ±29 

1866 ± 

2090 ± 59 

e 

0.345 ±0;03r 


0.265 ±0.015 

Stellar parameters 

Ra [R©] 

1 Qi 0 _|_ 0.095 
1.010 m 0.073 

1 OQG _(_ 0.066 
l.OUD It 0.073 

1.666 ±0.097 

Ma [Mq] 

1.300 ±0.060 

1 211 + 0127 

l.Zll It Q122 

1.350 ±0.120 

Pa [p©] 

n Kyn _|_ 0.120 
U.O/U m 0.092 

0-542 ±°:07g 

0.293 ±0.042 

logSA 

4 0-1 2 -U 0.061 
4.CJ1Z m 0.047 

4-287 ±0;0i 

4.126 ±0.042 

La [Lq] 

0.435 ±0.085 



Geometrical parameters 

a [au] 

0.037 ±0.001 

0.036 ±0.001 

0.037 ±0.001 


85.30 ±t^®° 

84.32 ±9:1? 

81.1±1.5 

b 

0.499 ±0:«4 

0.535 ±0:031 

0.752 ±0.133* 


*This parameter is not given in ISouthworthI 120121 ') but was 
calculated from Ra, a, i 


5 TRANSIT TIMING 

The mid-transit times were determined by applying the 
transit model created in Section |3] to the individual LCs 
using TAP. The initial parameters for the fitting were the 
ones given in Tabled The mid-transit time, as well as the 
flux slope and intercept were always a free parameter while 
the orbital period P and the eccentricity e were kept fixed. 
The best-model parameters and their uncertainties were set 
as Gaussian priors and i, u/Ra, Rh/RA, and the LD coeffi¬ 
cients of each transit were allowed to vary within this range. 
With this approach we assured that the best-model uncer¬ 
tainties are included in the error bars of the mid-transit 
time. Several LCs were observed in a filter that has not con¬ 
tributed to the template LC. In those cases the theoretical 
values ±0.5 was used as Gaussian prior for the LD coeffi¬ 
cients. Ten chains of a length of 10® steps were used for the 
MCMC analysis of each LC. Four transits were observed 
with more than one telescope. These LCs were fitted simul¬ 
taneously to increase the timing precision. The times have 
been converted into the barycentric Julian Date based on 
the barycentr ic dynamic time (BJ Dtdb) using the online 
converteiQ by lEastman et al.l 1 I 2 OIOII . 


5.1 Timing error estimates 

The observations at the University Observatory Jena and 
the YETI network yield 13 transit times. Looking at the 19 
individual LCs revealed some problems with several of the 
observations. 

® http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html 
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0.522 0.523 0.524 0.525 


Mid-Transit Time [JD-2455650] 

Figure 3. Distribution of the transit times obtained from 
the analysis of 185 LCs that were created as a transit time 
robustness check (see text for details) for the transit at epoch 
529 (2011 March 29). The thick solid line gives the result of 
TAP only for the initial LC, the dashed lines give the TAP 
error bars. The width of the transit time distribution, shown 
here as grey shaded area, is nicely reproduced by the TAP 
uncertainties (ratio between distribution width and TAP error 
bars = 0.8). 



2.555 2.56 2.565 2.57 2.575 2.58 

Mid-Transit Time [JD-2455630] 


Figure 4. Same as Fig. [3] but for the transit at epoch 521 
(2011 March 11). Given the smaller number of data points 
in the initial LC only 100 LCs could be created and anal¬ 
ysed. The TAP error bars clearly underestimate the width of 
the transit time distribution (ratio between distribution width 
and TAP error bars = 1.6). 


The transit from 2011 March 2 was observed with three dif¬ 
ferent telescopes but all of them only yielded a partial LC 
due to bad weather and technical problems. Since the egress 
is missing in all three cases, the shape of the transit cannot 
be determined properly. We excluded this data point from 
any further timing analysis. 

The data point from 2011 March 20 consist of three individ¬ 
ual transits that were modelled simultaneously and therefore 
has small error bars. But looking at the individual transits 
reveals that the transit time from the Calar Alto observa¬ 
tions differ by ~10min from the other two LCs. Since the 
reason for this difference stayed unclear (most likely the bad 
conditions like thin clouds, fog, and full moon during the ob¬ 
servation) we excluded the Calar Alto LC and modelled the 
remaining two simultaneously. 

The LCs from 2011 February 12 and 2011 March 11 suffer 


from outliers and bad weather conditions in the egress, re¬ 
spectively. In both cases we realized significant changes in 
the resulting transit times for different LC treatments like 
removing outliers or binning. Since these effects should be 
accounted for in the error budget we run a test to check the 
robustness of the transit times. From the original (raw) LC 
we first created a set of 100-200 LCs (depending on the ini¬ 
tial number of points in the LCs) by removing data points 
(every second, third, etc.), sigma clipping, using different 
binning factors, different artificial comparison stars, and re¬ 
moving trends. Then we modelled all LCs, determined the 
transit time and compiled a histogram. The width of the 
transit time distribution (represented by a Gauss function) 
gave the final timing error which were then compared to 
the outcome of tap. In most cases we could reproduce the 
uncertainties from the wavelet-based MCMC analysis (ratio 
between distribution width and tap error bars ~ 0.8-1.2) 
meaning the MCMC produces very robust error bars even 
for LCs with a very low quality. One example for a good 
agreement of the tap error bars and the distribution width 
is given in Fig. O Only for two LCs (2011 March 11 and 
2011 April 7) the tap error bars are by a factor of 1.5-2 
smaller than the distribution width as shown in Fig. [4] In 
both cases either the ingress or the egress is affected by the 
weather conditions which makes it difficult to recover the 
transit shape. 

As final uncertainties for the transit times, we always 
adopted the maximum of the tap error bars and the dis¬ 
tribution width. 


5.2 Transit ephemeris 

To exte nd the obse r vation al baseline, we included the data 
point of ijoshi et al.l (l2009t) to our timing analysis. Note, that 
this is not a single transit observation but the original pub¬ 
lished mid-transit time at epoch zero computed from many 
individual transits. Another transit observed with the Uni¬ 
versity of Haw aii 2.2 m (UH.2 . 2 m) t elescope on Manna Kea 
is published in I Johnson et al] (l2009l ). The mid-transit time 
for this observation was determined like for our measure¬ 
ments as described in Section [S] With these altogether 14 
mid-transit times (the partial transit at epoch 517 was ex¬ 
cluded) we re-calculated the transit ephemeris using an error 
weighted linear fit. The result is given in equation m, where 
E denotes the epoch (reduced = 0.64): 

rc[BjDTDB](-®) = (2454463.57688 -b E ■ 2.2437655) d , , 
±0.00047 ±0.0000010 


The orbital period P is 1.2 s longer and 10 times more precise 
than the one given in To shi et ahl (|2009|) but in agreement 
with ISouthworthl (|2012| ) and Blecic et al.l (l2013l ). 


5.3 Transit Timing Variations 

If we subtract the predicted transit times (calculated with 
the refined ephemeris) from the observed value, we obtain 
the 0-C. The results for the mid-transit times and the O- 
C are listed in Table m Fig. m shows the O-C diagram 
where the black solid line represents the refined ephemeris. 
To search for a periodicity in the O-C diagram we com¬ 
puted the generalized Lomb-Scargle periodogram (gls; 
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IZechmeister fc Kursterir2009l ~). The periodogram is shown in 
Fig. |6] where the highest peak (Pttv = 26.33 ± 0.09 epochs, 
power of 0.66) corresponds to a False-Alarm probability 
(FAP) of 44.7%. Within our data set we could not detect 
any evidence for TTV. 

Residuals in transit timing allow us to place constraints on 
the properties of any hypothetical additional planet in the 
system. A set of synthetic O-C diagrams for WASP-14 b in 
the presence of a fictitious per turbing planet w as generated 
with the MERCURY 6 package (IChamber j|l999l l and the im¬ 
plemented Bulirsch-Stoer integrator. The mass of this sec¬ 
ond planet was set at 0.5, 1, 5, 10, 50, 100, and 500 Msarth 
(Earth masses). The initial semimajor axis was iterated from 
0.01 to 0.12 an with a step of 2 x 10“® an. The system was as¬ 
sumed to be coplanar. For WASP-14 b, the initial argument 
of periastron, ujh, and orbital eccentricity, es were taken from 
Blecic et al. (2013), and the mean anomaly was set at 0°. 
For the fictitious planet, its argument of periastron was set 
at the value equal to ojh at the beginning of each simula¬ 
tion. Calculations were done for three values of the orbital 
eccentricity of the fictitious planet, Cp, set to 0.0, 0.1, and 
0.2. The initial mean anomaly was shifted by 180°. The inte¬ 
gration for each planetary configuration covered 2250 d (i.e. 
1000 orbit of WASP-14 b - the time span of transit timing 
observations). We calculated the rms of the signal at Pttv 
as rmsTTV = Amax/v^, where Amax is an amplitude of this 
signal. We derived rmsTTV = 37 s. Analogously, rms of the 
residuals from a linear ephemeris was calculated for each set 
of simulated observations. Then, for each orbital distance, 
we determined the range of planet masses in which the value 
of rmsTTV fell. An upper mass of the fictitious planet at the 
detection limit was found by linear interpolation for masses 
below 500 MEarth- If the TTV sinal was found to be gener¬ 
ated by a more massive body, the limiting mass was extrap¬ 
olated using a linear trend as fitted to 100 and 500 MEarth- 
Exemplary results for Cp = 0.1 are illustrated in Fig. [7] The 
timing technique allows us to probe the Earth-mass regime 
close to MMRs. Most of orbits located between the inner 
1:2 and outer 2:1 orbital period commensurabilities, were 
found to be highly unstable and planetary close encoun¬ 
ters or planet ejections occurred during the relatively short 
time of integration. Thus, planetary conhgurations which 
are placed in this space are highly unlikely to exist. 


25 
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Epoch 

Figure 5. The O—C-diagram of WASP-14 b. The black filled 
and open (with dashed error bars) symbols denote the complete 
and the partial transits, respectively. Note that the transit at 
epoch 669 (2012 February 6) is considered here as ‘complete’ 
since the ingress as well as the egress, a nd therefore the tr ansit 
shape, are intact. T he data points from Ijoshi et alj and 

Ijohnson et alJ l|20Q9h are shown in grey. The solid line represents 
the updated ephemeris given in equation ®. 
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Figure 6. Generalized Lomb-Scargle periodogram (top panel) 
and window function (bottom panel) for the O—C dia¬ 
gram of WASP-14 b. The highest peak with a period of 
Pttv = 26.33 ± 0.09 epochs at a power of 0.66 shows a FAP of 
44.7%. 


6 CONCLUSIONS 

WASP-14 b is, for several reasons, a very interesting target. 
First, because of its very high mass it is one of the densest 
exoplanet with an orbital period shorter than three days. 
Furthermore, despite of its close-in orbit WASP-14 b has a 
rather high eccentricity. Interestingly, there is a strong ten¬ 
dency that massive planets on close-in orbits have ellipti¬ 
cal orbits. Approximately 58% (7 out of 12, exoplanet.eu, 
2014 November 12) of transiting exoplanets with masses 
greater than M = 5 Mjup and periods less than 10 d have 
e 7 ^ 0 while only 20% (53 out of 265, exoplanet.eu, 2014 
November 12) of the less massive planets show a signifi¬ 
cant eccentricity. This may indicate that there are two dis¬ 
tinct types of exoplanets. Another hint is the distribution in 
the period-mass diagram for close-in exoplanets (P < 10 d 
and M < 15Mjup) which is shown in Fig. [S] Two dis¬ 


tinct areas are clearly identified. Most of the planets have 
a mass below 5Mjup. So far, no planets were discovered in 
the mass range between 4.5 and 6.5 Mjup (The only point 
in this area belongs just to an upper mass limit of WASP- 
33 b, M < 4.59 Mjup), while several objects were found with 
masses ^ 6.5 Mjup, which are often referred to as hot super- 
Jupiter. This clear distinction could indicate a physical dif¬ 
ference in the planet formation and evolution. Since only a 
handful of hot super-Jupiter known so far, this massive area 
is insufficiently explored to make general statements. 

We observed 19 LCs of WASP-14 b with six telescopes at 
five different observatories within the TTVQYETI collabo¬ 
ration. Due to the weather conditions and observations done 
with small telescopes the LCs are of highly variable quality. 
All transits are shown in Figsl^and llOl grouped according to 
their quality to LC with rms > 4 mmag and rms < 4 mmag. 
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Table 7. Transit times for all observed transits of WASP-14 b including the publicly available transits. The O—C was calculated with 
the ephemeris given in equation 


Date 

Epoch 

Tc (BID) 

O-C (min) 

Ref. 


0 

2454463.57657 ±0.00053 

-0.44 ±0.76 

Joshi et al. ('2009') 

2009 Apr. 1 

205 

2454923.54564 ±0;00310 

-4.55±|:|0 

This work 

2009 Apr. 19 

213 

2454941.50043 ±g;gg||| 

2-16 ±?:i 

This work 

2009 May 7 

221 

2454959.45027 ±g;gj?gg 

1 7 c -L 14.40 
10.94 

This work 


223“ 

2454963.93776 ±0;00°®g 

1 09 + 0.99 
i.uy in 

Johnson et al. f2009) 

2011 Feb. 12 

509 

2455605.65997 ±g;g)J|fg 

4-83 ±111 

This work 

(2011 Mar. 2 

517'’ 

2455623.59980 

-5.53 ±y 6 

This work)* 

2011 Mar. 11 

521 

2455632.57247 ±0:00306 

-3.01 ±1|0 

This work 

2011 Mar. 20 

525'’ 

2455641.54831 ±0;00055 

- 0 . 86 ± 0 ;l| 

This work 

2011 Mar. 29 

529'’ 

2455650.52899 ±g;ggg|| 

n 24 -h 0-84 

m Q g4 

This work 

2011 Apr. 7 

533 

2455659.50506 ±g gg||| 

1.69 ±1|0 

This work 

2011 Apr. 16 

537 

2455668.47584 ±0-00503 

-4.48 ±111 

This work 

2012 Feb. 6 

669'’ 

2455964.65659 ±g;gg(|g 

0.86 ±102 

This work 

2012 Feb. 24 

677 

2455982.60621 ±0-00090 

0.13 ±li 

This work 

2013 Apr. 30 

869 

2456413.40914 ±0-00163 

0-07 ±i;i 

This work 


*Data point was excluded from timing analysis. 

“Transit was re-analysed with TAP. 

^The LCs have been combined during the MCMC analysis to improve timing precision. 



Pp/Pb 


Figure 7. Upper mass limit for a fictitious additional planet in 
the WASP-14 system, based on the timing data set, as a func¬ 
tion of the orbital period of that planet, Pp. The greyed area 
shows unexplored configurations because they are below a detec¬ 
tion threshold of the timing data set. 



0 2 4 6 8 10 


orbital period [d] 

Figure 8 . The period-mass diagram for close-in exoplanets (P < 
10d and M < 15Mjup). Only 12 planets are more massive than 
6.5Mjup. No planets were found in he region between 4.5 and 
6.5Mjup (except WASP-33 b which is just an upper mass limit 
M < 4.59 Afjup, upper limits are marked with an arrow) 


respectively. 

From the simultaneous LC modelling of our four best LCs 
we could determine planetary, stellar and geometrical prop¬ 
erties of the system. Our v alues are in agreement with the 
values in ijoshi et al.l ll2009l l but differ sign i ficant ly from the 
physical properties given in ISouthworthI (l2012l l. This dis¬ 
crepancy, however, is not physical since it is only caused by 
the interpretation of the final results bv ISouthworthI ||2012| . 
averaging of two different sets of best-fitting parameters for 
two individual LCs). Since the two high precision transit LCs 
available in literature were found to be inconsistent with 
each other our findings are from great importance for the 
determination of the system parameters. 

Including the two publicly available data points, altogether 


14 mid-transit times (the partial transit at epoch 517 was 
excluded) were used in the transit timing analysis. To inves¬ 
tigate the error budget we ran a test to check the robustness 
of the transit times. We found, that if the ingress or egress, 
and hence the shape of the transit, is missing or show outlier 
and/or large scatter even the uncertainties determined with 
a wavelet-based MCMC may be underestimated. 

We found no significant periodic signal in the O-C diagram. 
The strongest period at ~26 epochs has a FAP of 44.7%. 
Hence, there is no evidence for TTV in the system. From 
our three-body-simulations, we can exclude even earth-mass 
perturbers in some resonance-orbital configurations. 
Measurements of the Rossiter-McLaughlin effect found a 
misalignment between the stellar spin axis and the orbital 
axis of the planet which cannot be explained by planet-disc 
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migration theories ie.g. iGaudi fc WinrJl2007l '). Another pos¬ 
sible mechanism that may be responsible for such massive 
clos e-in planets is gravitational scatt ering by larger plan¬ 
ets dWeidenschilling fc Marzaril Il996l l. This planet-planet 
scattering scenario is one way to explain the properties of 
WASP-14 b. 

The significant eccentricity indicate either that the tides 
did not had sufficient time to influence the planetary orbit, 
which would support the planet-planet scattering theory or 
that the tidal effects in planetary systems are weaker than 
expected. Long-term follow-up studies of WASP-14 will add 
stricter constraints on these theories. 
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Transit Timing of WASP-lf b 
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Figure 9. LCs of WASP-14 b with an rms > 4 mmag. The date of observation, observatory, filter, pnr, and the rms of the fit are indicated 
in each individual panel 
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Figure 10. The same as Fig. [9| but for the higher quality LCs with an rms<4mmag. The four best quality LCs (in terms of pnr) that 
were used to create the template are marked with an asterisk. 
















































